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Abstract 

We study a set of exactly soluble net spin models. There exist two kinds 

of ground state, one is a complete dimerized state, and the other one is the 

ground state of corresponding spin-1 model. For the excitation gap, various 

phases were discovered and determined. 
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It is well known that a spin-1/2 antiferromagnetic chain with nearest-neighbor coupling 
Ji and next neighbor coupling J2 could dimerize in the presence of frustration, as exemplified 
by the exactly soluble one- dimensional Majumdar-Ghosh model where the exact (twofold 
degenerate) ground state is a simple product of singlet dimers. The elementary excitation 
can be constructed as a pair of unbound spins above the completely dimerized state 0. 
This model also can be considered as a two-chain with one diagonal or a zigzag spin ladder. 
Analytical and numerical studies f^-^l of the model show a transition from a gapless phase 
when J2 < J2C to a gapped phase when J2 > J2c- Recently, the Majumdar-Ghosh model 
was studied by the density matrix renormalization group (DMRG) approach. The DMRG 
study of Chitra et al. considered a more general model which included a dimerization 
term. White and Affleck studied the zigzag spin ladder by using the bosonization and the 
DMRG methods. In particular, they pointed out a potential important connection between 
the zigzag chain and the generalized Kondo lattice. In addition, a family of spin-ladder 
models were studied by Kolezhuk and Mikeska p] that exhibit non-Haldane spin-liquid 
properties as predicted by Nersesyan and Tsvelik 

In this Letter we consider a set of spin-1/2 antiferromagnetic Heisenberg models defined 
on a double layer as shown in Fig. 1. Each layer has M x L sites and it connects to the 
other one not only perpendicularly, but also diagonally. L is a measure of the dimensionality 
which could be any integer varying from 1 (the net spin ladder) to M (the net spin layers). 
Throughout this work, we use S to represent spins on the lower layer, and S' on the upper 
layer. Open boundary conditions are imposed. 

Under certain combination of coupling constants, we show that for general L and M 
the model could be solved exactly. First, the model have two kinds of ground state, one 
is a completely dimerized state, and the other one is the ground state of the spin-1 model 
defined on the single layer. Second, the model exhibits rich excitation phases, depending 
on coupling constants and dimensionality L. As we shall show below, these properties are 
quite different from the existing models such as the Majumdar-Ghosh model 

The model Hamiltonian is: 
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(1) 



fc,Z=l 



with Ji, J2, J3, J4 > 0, and |S| = |S'| = S = 1/2. Both Ji and J2 favor local antiferromag- 
netic ordering while J3 and J4 represent frustration effects. 

For any integers M and L, we found that a complete dimerized state: 



where [i,j] denotes the normalized singlet -^{aiPj — Piaj) with a and P representing the 
usual up and down single spin eigenf unctions, respectively, is an eigenstate of Eq. (||) when 
we impose the following constraint: 



To prove our statement, we first show that this is true for the 4-spin plaque case, 
{Si^i, S'l/^i', 82,1', S'2',1'}, simplied as {Si, S'l', S2, S'2'}, which is indicated in Fig. 1. The 
model Hamiltonian commutes with total spin Stot so one can use it to characterize eigen- 
states. In the sub-space of zero 2;-component magnetization there are six eigenstates and 
among these six eigenstates there are two singlets which are determined by the following 
equations: 



[1, l'][2, 2'] ■ ■ ■ [M, M'] ■ ■ ■ [iV, N'], N = ML , 



(2) 



2J2 — t/3 + t/4 • 



(3) 



HlP^ = ^(2J2 - J3 - JMd + (Jl - 2J2 - J3 - J4)V^2 



(4) 



where 




(5) 
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is the completely dimerized state it becomes an eigenstate of H with eigenvalue — 3Ji and 
another singlet state ip2 with eigenvalue E2 = Ji — 4 J2 when Eq. (^) holds. 

Also, we obtained other four eigenvalues according to their total spin Stot, among them 
there is a triplet with eigenvalue: 

E3 = -iJl + J2) + liJ3 + J4) ■ (6) 

We know that the formation of dimers [1, 1'] and [2, 2'] is due to antiferromagnetic cou- 
pling Jl. Without antiferromagnetic couplings J3 and J4, quantum fluctuation due to J2 
will kill both dimers. As long as condition Eq. @ holds, the effect of J2 on the dimers will 
be canceled exactly by J3 and J4. This is also true no matter how many plaques we put 
together and thus ipD is the eigenstate of Hamiltonian (P with eigenvalue Ed = — |MJi. 

It is also obvious that the dimerized state ipn may not be the ground state for general 
coupling parameters {^1,^2,^3,^4} due to the fact that the antiferromagnetic couplings 
J2, J3, Ji which connect one plaque to another also introduce interactions between plaques. 
Since in the limit of Ji 00, ip^ must be the ground state, and in the limit of Ji 0, ifju 
cannot be the ground state, therefore we expect that there exists a critical value of Jic such 
that when Ji > Jic the ground state is completely dimerized. 

We should mention that for the case of L = 1, the spin net model (we call it net spin 
ladder) goes back to the Majumdar-Ghosh model ^ when J2 = 0.5Ji, J3 = Ji, and J4 = 0. 
Several studies on the generalization of the M-G model, such as the zigzag model , were 
carried out before. 

To proceed, we study the case of L = 1 with coupling constants J3 = J4 = J2 in details 



as an illustration. This special case has been recently studied by several groups |[T0|- [15|] . We 
use it as an example to present our approach which leads us to solve the whole class of net 
spin models. 

We can rewrite the Hamiltonian as 

Q M-1 M-1 

H = -2MJi- + J,Y,al + 2J2 '^fc ■ "^k+i , (7) 

^ k=l k=l 



where crk = Sk + S'k. This Hamiltonian describes a chain of M spins with either \crk\ = 
(singlet) or 1 (triplet). 

With this form of H, it is easy to see that ipD is an eigenstate of H by using the relation 
(Sfc + S'k)[k, k'] = 0, hence 



/ M-l M-l \ 

Ji ^ + 2J2 J2 ■ o-k+i li^D >= 

\ k=l k=l J 



(8) 



and the eigenvalue is E^, = — |MJi. 



To lower energy, the squared term, Jicr^, favors to be in the |crfc| = state, while the 
exchange term, J2CA: ■ cT/t+i, favors to be in the |crfc| = 1 state. If all pairs + S'^ are in 
the (7 = 1 state then the model is equivalent to the spin-1 chain with coupling constant 2J2, 
apart from a constant One can see from Eq. (0) that the eigenstates of the spin-1 chain 
could also be the eigenstate of the net spin ladder model. In fact, as Ji decreases further, 
the ground state of the net spin ladder model becomes the ground state of the spin-1 chain. 

Moreover, we can use this criterion to determine the transition point Jic. When Ji = Jic, 
< H >= Ed = — f MJi, where the expectation value is taken in the ground state. We obtain 

M-l M-l 

< Jic J2 ^1 + 2'^2 J2 ■ CTk+l >= , (9) 

k=l k=l 

i.e., 2JicM = ground state energy of M-site spin-1 chain with coupling constant 2J2. Using 



the existing estimates ||16[, we obtain 

Jic = I.4OI5J2 • (10) 

To confirm our analysis, we determine the critical value Jic by using the exact diagonalization 
technique for several clusters. Our numerical answer is Jic = 1.402 ± 0.001, the same 



as the above within numerical error. 

We next study excitation spectrum of the model. One can create a single magnon simply 
by breaking a dimer: -^{uiPj — Pittj) — > aiaj (or PiPj or -^(aiPj + Piaj).) It can be shown 
that state 

V^™ = [1,1'][2,2']---(T,T)---[M,M'] , (11) 



is also an eigenstate of the net spin ladder model with total spin Stot = 1- The energy gap 
in this case is exactly A^j = 2 Ji and the transition is from singlet to triplet. This excitation 
energy is independent of other antiferromagnetic couplings as long as the condition holds 
and thus it is possible that there exists states with lower excitation energy. 

Analyzing results for the single plaque, we found that when Ji < 2 J2, E2 = Ji — 4J2 is 
lower than = — Ji of Eq. (P). The former is in the singlet {Stot = 0) sector. Interesting 
enough, this critical value (Ji/J2)c = 2 is independent of lattice size. This could be easily 
understood by using the Hamiltonian of Eq. (|^ . Such excitation corresponds to the situation 
where two adjunct dimers are broken (so they are in the spin 5 = 1 state) while the rest 
of dimers remain unbroken, regardless the chain length. The two spin 1 sites form a singlet 
with energy — 4J2 so the excitation energy, which corresponds to the transition from the 
dimer singlet to another singlet, is A^^ = 4Ji — 4J2. A^^ < A^^ when Ji < 2J2. 

What will happen as one decreases Ji further? We found that the first excited state 
remains to be the one we just described for the whole region of 2 J2 > Ji > Jic- Moreover, 
this state is (M — l)-fold degenerate. States consisting of longer segment of broken dimers 
are high in energy. 

Figure 2a summarizes our findings. Here gap A is defined as the difference between the 
first excited state energy and the ground state energy, E{1) — E{0). There are three regions, 
depending on parameter J2/J1: 

Region I , Ji > 2J2, A/Ji = 2, independent of J2 and it is singlet to triplet. 

Region II , 2J2 > Ji > Jic, A/Ji = 4(1 — J2/J1) decreases linearly from 2 to 4(1 — J2/ Jic) = 
1.1459 at Ji = Jic+. The excitation is from singlet (dimers) to another singlet. In 
both regions I and II, the completely dimerized state is the ground state. 

However, right at Ji = Jic, the completely dimerized state is degenerate with the 
ground state of spin-1 chain. 

Region III , Ji > Jic, it is the Haldane gap phase. At Ji = Jic—, A/Ji = 0.5851. 
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We now turn to the double layer case {L = M and = LM). Similarly, we get 

q Af N 

H = -2NJ,- + J,Y,(tI + 2J2 "^k-cTi, (12) 

k=l <k,l> 

where cr^ = S^,; + S'^,; and < k,l > refer to nearest neighbors in two dimensions. Again, 
when all |cr| = 1, this is nothing but the two-dimensional antiferromagnetic Heisenberg 
model with coupling constant 2J2, apart from a constant. 

With the same approach we had described above, it is easy to see that a complete 
dimerized state where all cr are in singlet state is the eigenstate of the model and it is the 
ground state when Ji is sufficiently large. The critical point at which the complete dimerized 
state becomes the excited state is determined by, 

= -1^ (^k -o-i . (13) 

<k,l> 



Using the existing estimates |T8[, we obtain 

Jic = 2.3323J2 . (14) 
Two interesting observations are right in hand: 

(1) The existence of the critical point is similar to that for the double layer layer antifer- 
romagnetic Heisenberg model without diagonal coupling (J3 and J4 set to be zero) terms. 
There has been a lot of studies on that model, by the spin-wave theory, perturbation theory, 
bond-operator approach, etc. (Ref. [19] and references therein.) The transition point at 
which gap vanishes is J1/J2 = 2.54, close to the value of Jic we got here, although with 
different model. 

(2) We can show that it is still true when J1/J2 < 2, the first excitation is a triplet. 
Moreover, after J1/J2 < 2.3323, the system is gapless because the existence of long-range- 
order in the two-dimensional spin-1 Heisenberg model, as shown rigorously long time ago 

20| , |2T[] . Therefore, there exists no transition from the dimer singlet to another singlet as 



what we have seen for the net spin ladder case. Note that the numerical accuracy of 2.3323 
is unimportant here, as long as the ground state energy per site is lower than — 2J2, (which 



must be true because — 2J2 just corresponds to the Neel state!), our assertion holds. We 
plot the excitation gap as function of J2/J1 in Fig. 2b. 

What will happen for cases of 1 < L < 00? We have studied two cases, L = 2 and L = 3. 
For L — 2, the critical value Jic is determined by the corresponding spin-1 model defined on 
the 2 X M lattice. Table I. hsts the ground state energy per site as well as per bond. Total 
number of lattice sites is A'^ = 2 x M and total number of bonds is A*";, = 3 x M — 2 when 
open boundary is used. 

It is clear that, in the limit of A?" — > 00, we have 

E/NJ2 = {E/N,){N,/N) = {3/2){E/N,) (15) 
> (1.5) * (-1.3) = -1.95 > -2.0. 

Thus, Jic < 2J2 and the excitation spectrum for the case of L = 2 is similar to that of 
L= 1. 

Similarly, we list the ground state energy per site as well as per bond in Table II for the 
case of L = 3. 

Because we do not have enough data due to lack of computer power, we cannot confidently 
establish a lower bound for the ground state energy in the thermodynamics limit. Although 
a second order perturbation theory calculation gives 

617 630 
= -315 > 316 = • 

it is too close to be definitive. On the other hand, the second order perturbation theory 
calculation for the L — 2 case gives 

E/NJ2 = -| > -2 , (17) 
o 

which is more definitive. 

Therefore, we cannot rule out a possibihty of dimensional crossover from 1-D to 2-D 
occurs at L = 3. It is also equal probable that Jic < 2J2 but the difference between Jic and 
2J2 becomes exponentially small when L increases from 3. 



For other combinations of couplings, one can only determine the critical points and 
eigenvalues of the model numerically. However, qualitative behaviors should be the same as 
what we have shown here. Note that the energy gap 2Ji is quite large (critical values of J2 
are smaller than Ji) so the phase diagram should be similar to Fig. 2. Also we note that 
one obtains a critical value of J1/J2 = 2.54 when J3 = J4 = is very close to that when 
J3 = J4 = J2 = (1/2.33) Ji. This shows that our results are insensitive to the values of J3 
and J4. Essential physics is determined by the two couplings: J2 and Ji. 

In summary, we have studied a class of net spin models, defined by the Hamiltonian 
Eq. (|l]), and controlled by three antiferromagnetic couplings: J2/ Ji, J3/J1, J4/J1, and 
dimensionality L. When 2J2 = J3 + J4, the completely dimerized state ipD is an eigenstate 
of these models and it is also the ground state if Ji is sufficiently large. To fully understand 
the model, we have studied the case of J2 = J3 = -h extensively and obtained phase diagrams 
for both 1-D (ladder, L = 1) and 2-D (double layer, L = M) cases. For the ladder case at a 
particular point in the parameter space: J2 = jJi, J3 = Ji, and J4 = this model returns 
to the Majumdar- Ghosh model But its behavior is quite different from that of the 
Majumdar- Ghosh model. We also show qualitative different behavior between the ladder 
and the double layer cases. Moreover, we have studied the cases of L = 2 and L = 3 and 
discussed the possibility of crossover from 1-D to 2-D. In addition, we made connection with 
other studies of double layer system. 
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FIGURES 




FIG. 1. The net spin layers of 2ML spins, where thick solid lines represent singlet dimers. 
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Fig. 2a 
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Fig. 2b 
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FIG. 2. Excitation gap as function of coupling parameters for the net spin ladder (a), and 
double layer (b) models. 
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Table I 



Lattice 


E{0)/N 


E{0)/Nb 


2x2 


-1.50000000 


-1.50000000 


2x3 


-1.61633441 


-1.38542949 


2x4 


-1.68229677 


-1.34583741 


2x5 


-1.72097660 


-1.32382816 


2x6 


-1.74706131 


-1.31029598 


2x7 


-1.76571007 


-1.30104952 


2x8 


-1.77973127 


-1.29435001 


Table II 


Lattice 


E{0)/N 


E{0)/N, 


3x2 


-1.61633441 


-1.38542949 


3x3 


-1.71359964 


-1.28519973 


3x4 


-1.80819446 


-1.27637256 


3x5 


-1.83887667 


-1.25377955 



